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Abstract 

We present an exact and Monte Carlo renormalization group (MCRG) study of semiflexible 
polymer chains on an infinite family of the plane-filling (PF) fractals. The fractals are compact, 
that is, their fractal dimension df is equal to 2 for all members of the fractal family enumerated 
by the odd integer b (3 < b < oo). For various values of stiffness parameter s of the chain, on the 
PF fractals (for 3 < b < 9) we calculate exactly the critical exponents v (associated with the mean 
squared end-to-end distances of polymer chain) and 7 (associated with the total number of different 
polymer chains). In addition, we calculate v and 7 through the MCRG approach for b up to 201. 
Our results show that, for each particular b, critical exponents are stiffness dependent functions, 
in such a way that the stiffer polymer chains (with smaller values of s) display enlarged values 
of v, and diminished values of 7. On the other hand, for any specific s, the critical exponent v 
monotonically decreases, whereas the critical exponent 7 monotonically increases, with the scaling 
parameter b. We reflect on a possible relevance of the criticality of semiflexible polymer chains on 
the PF family of fractals to the same problem on the regular Euclidean lattices. 
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I. INTRODUCTION 



The self-avoiding walk (SAW) is a random walk that must not contain self-intersections. 
It has been extensively studied as a challenging problem in statistical physics, and, in par- 
ticular, as a satisfactory model of a linear polymer chain The pure SAW is a good model 
for perfectly flexible polymer, where we ignore the apparent rigidity of real polymer, and, 
consequently, to each step of SAW we associate the same weight factor (fugacity) x. In most 
real cases, the polymers are semiflexible with the various degree of stiffness. To take into 
account this property of polymers, in the continuous space models the stiffness of the SAW 
path is modeled by constraining the angle between the consecutive bonds of polymer, while 
in the lattice models, an energy barrier for each bend of the SAW is introduced. The lattice 
semiflexible SAW model (also known as persistent or biased SAW model), has been studied 
some time ago in a series of papers [2j, with a focus on the so-called rod-to-coil crossover. 
Afterwards, it was modified in various ways, in order to describe relevant aspects of different 
phenomena, such as protein folding jjj, J], adsorption of semiflexible homopolymers j^], tran- 
sition between the disordered globule and the crystalline polymer phase {(], 7], behavior of 



semiflexible po 

systems 
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ymers in confined spaces [si, 



or influence of an external force on polymer 



In spite of numerous studies, a scanty collection of exact results for semiflexible polymers 
has been achieved so far, even for the simplest lattice models. A few cases in which some 
properties of semiflexible SAW can be studied exactly are: directed semiflexible SAWs on 



regular lattices 
fractal lattices 



14 L and semiflexible SAWs (with no constraints on the direction) on some 
161 ] . In particular, exact values of the end-to-end critical exponent v and 



the entropic exponent 7 were obtained for these models, and it turned out that in some cases 
critical exponents are universal, whereas in other cases they depend on the stiffness of the 
polymer chain. Universality arguments, as well as results of approximate and extrapolation 
methods for similar models suggest that critical exponents on regular (Euclidean) lattices 
should not be affected by the value of the polymer stiffness. On the other hand, it is not 
known what are the effects of rigidity on the critical behavior of SAWs in nonhomogeneous 
environment. In order to explore further this issue, in this paper we perform the relevant 
study on the infinite family of the plane-filling (PF) fractal lattices [l?], Q], which allow 



for an exact treatment of the problem. These fractals appear to be compact, that is, their 
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oo characteristics of these fractals approach, via 
Q, properties of the regular 2D lattice. 
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221 ]. as well as 



261 ]. we calculate critical 



fractal dimension df is equal to 2. Members of the family can be enumerated by an odd 
integer b (3 < b < oo), and as b 
the so-called fractal-to-Euclidean crossover 

By applying the exact real-space renormalization group (RG) method 
Monte Carlo renormalization group (MCRG) method [2C 
exponents v and 7. We have performed our calculations for as many as possible members 
of the fractal family, for various degree of polymer stiffness, in order to study consequent 
stiffness dependence of the critical exponents, as well as to see the asymptotic behavior of 
the exponents in the fractal-to-Euclidean crossover region. 

This paper is organized as follows. We define the PF family of fractals in Sec. [TT1 where 
we also present the framework of our exact and MCRG approach to the evaluation of the 
critical exponents v and 7 of stiff polymers on the PF fractals, together with the specific 
results. In Sec. IIHI we analyze the obtained data for the critical exponents, and present an 
overall discussion and pertinent conclusions. 



II. SEMIFLEXIBLE POLYMERS ON THE PLANE-FILLING FRACTAL LAT- 
TICES 

In this section we are going to apply the exact RG and the MCRG method to calculate 
asymptotic properties of semiflexible polymer chains on the PF fractal lattices. Each member 
of the PF fractal family is labelled by an odd integer b (3 < b < 00), and can be constructed 
in stages. At the initial stage (r = 1) the lattices are represented by the corresponding 
generators (see Fig. [Tj). The rth stage fractal structure can be obtained iteratively in a 
self-similar way, that is, by enlarging the generator by a factor b r ~ 1 and by replacing each 
of its segments with the (r — l)th stage structure, so that the complete fractal is obtained 
in the limit r — > 00. The shape of the fractal generators and the way the fractals are 
constructed imply that each member of the family has the fractal dimension df equal to 
2. Thus, the PF fractals appear to be compact objects embedded in the two-dimensional 
Euclidean space, that is, they resemble square lattices with various degrees of inhomogeneity 
distributed self-similarly. 

In order to describe stiffness of the polymer chain, we introduce the Boltzmann factor 
s = e~ A / kBT , where A is an energy barrier associated with each bend of the SAW path, and 
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(a) (b) 

FIG. 1: (Color online) (a) The first three fractal generators (r = 1) of the plane-filling (PF) family 
of fractals, (b) The fractal structure of the 6 = 5 PF fractal at the second stage of construction, 
with an example of a piece of a possible SAW path (thick line). The full dots represent the turn 
points of the walker (that is, the bends of the SAW path), to which we associate the Boltzmann 
factor s = e~ A / fesT , where A > is the energy of SAW bend. Thus, for example, the presented 
SAW configuration should contribute the weight x 97 s 62 in the corresponding RG equation (more 
specifically, in equation (12. 6ft for r = 0). 

ks is the Boltzmann constant. For 0<s<l(0<A< oo) we deal with the semifiexible 
polymer chain, whereas in the limits s — 1 (A = 0) and s = (A = oo) the polymer is a 
flexible chain or a rigid rod, respectively. If we assign the weight x to each step of the SAW, 
then the weight of a walk having A" steps, with Nb bends, is x N s Nb , and consequently, the 
general form of the SAW generating function can be written as 

oo N-l 

G(x,s) = E C(N,N b )x N s N ", (2.1) 

N=l N b =0 

where C(N Nb) is the number of Af-step SAWs having Nb bends. For large N, it is generally 
expected [l| that the total number C(N,s) = En^2 C(N, N b )s Nb of N-step SAW displays 
the following power law 

C(N, s) ~ /^AH- 1 , (2.2) 

where 7 is the entropic critical exponent, and \x is the connectivity constant. Accordingly, 
at the critical fugacity x c = l//x(s), we expect the following singular behavior of the above 
generating function 

G sing ~ (x c - xY 1 . (2.3) 



On the other hand, due to the self-similarity of the underlying structure, an arbitrary 
SAW configuration on the PF fractals can be described, by using the three restricted gen- 
erating functions A^ r \ 5" and (see Fig. [2]), which represent partial sums of statistical 
weights of all feasible walks within the rth stage fractal structure for the three possible kinds 
of SAWs. One may verify that, for arbitrary b, the generating function G(x, s) is of the form 



00 1 f r 
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g 2 (Bl*- 1 \8)\C<r- 1 '>] , (2.4) 



where the coefficients gi(B^ r ~ l \ s) and g 2 (B^ r ~ 1 \ s) are polynomials in _g( r_1 ) and s. This 
structure for G(x, s) stems from the fact that all possible open SAW paths can be made 
in only two ways, using the rth order structures. The functions A^'\ B^ and appear 
to be parameters in the corresponding recursion (renormalization group) equations, which 
have the form 

A^ = a{B { - r \s)A^ , (2.5) 
EH- 1 ) = 6(5« )S ) , (2.6) 
C (r+1) = Cl (B {r \s)(A {r) f + c 2 (B {r \s)C {r \ (2.7) 

where the coefficients a(B^ r \ s), b(B^ r \ s) and Ci(B^ r \ s) (i = 1, 2), are polynomials in terms 
of and s, and do not depend on r. The established RG transformation should be 
supplemented with the initial conditions: A^ = ^Jx , B^> = x, and = 0, that are 
pertinent to the fractal unit segment. 

The basic asymptotic properties of SAWs are characterized by two critical exponents v 
and 7. The critical exponent v is associated with the scaling law (R%) ~ N 2u for the mean 
squared end-to-end distance for A-step SAW, while the critical exponent 7 is associated 
with the total number of distinct SAWs described by f!2.2p . for very large N We start 
by applying the above RG approach to find the critical exponent v for semiflexible polymers 
on PF fractals. We shall first present the corresponding exact calculation, and then we shall 
expound on the MCRG approach. To this end, we need to analyze (12 .6p at the corresponding 
fixed point. It can be shown that b(B^ r \s) in (J2J| is a polynomial 

B' = J2 p ( N > N b )B N s Nb , (2.8) 

N,N b 




A 



(r) 





FIG. 2: (Color online) A diagrammatic representation of the three restricted partition functions 
for an rth order of the fractal construction of a member of the PF family. The fractal interior 
structure is not shown. Thus, for example, represents the SAW path that starts somewhere 
within the fractal structure and leaves it at its upper right link to rest of fractal. 



where we have used the prime for the rth order partition function and no indices for the 
(r — l)th order partition function. Knowing the RG equation (12.81) . value of the critical 



exponent v follows from the formula 



22] 



ln6 
In \ v 



(2.9) 



where X u is the relevant eigenvalue of the RG equation (12. 8p at the nontrivial fixed point 
< B*(s) < 1, that is 



A,. 



dB' 
~dB 



B* 



(2.10) 



Consequently, evaluation of v starts with determining the coefficients P(N, N^) of (12.81) and 
finding the pertinent fixed point value B*(s), which is, according to the initial condition 
B^ = x, equal to the critical fugacity x c = 

We have been able to find exact values of P(N, iV&) for 3 < b < 9. For the first two 
fractals (6 = 3 and 5) the RG equation (I2.8P has the form 



B' = B 3 + 2B 5 s 



5„4 



(2.11) 



B' = B 5 + 125V + 2EV + 12EV + QB 9 s 8 

AB n s 6 + 8B u s s + 2B 13 s s + 4B 13 s w + 2B 15 s 12 , 

(2.12) 

respects, while for 6 = 7 and 9 they are disposed within the Electronic Phys.cs Auxiliary 
Publication Service (EPAPS) [27|. Knowing P(N,N b ), for a given b, we use Eqs. (12. 8ft - 
(I2.10p to learn B*(s) and the critical exponent v{s). For the 6 = 3 PF fractal the critical 
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fugacity B* and the critical exponent v can be obtained in the closed forms as functions of 
the stiffness parameter 



B U (s) ^™-\ (2.13) 

*-<•> = i^ragnj • (2.i4) 

while for b = 5, 7 and 9 they can only be calculated numerically. We have chosen the set 
of six values for the polymer stiffness parameter (s = 1,0.9,0.7,0.5,0.3 and s = 0.1) and 
the obtained exact values are presented in the Tables [I] and [III (together with the results 
obtained by the MCRG method). 

To overcome the computational problem of learning exact values of P(N, Nb), for fractals 



with b > 11, we apply the Monte Carlo renormalization group method (MCRG) [261 ] . The 
essence of the MCRG method consists of treating B', given by (12. 8p . as the grand canonical 
partition function that accounts for all possible SAWs that traverse the fractal generator at 
two fixed apexes. In this spirit, (12.81) allows us to write the following relation 

dB' B' 



dB B 

where (N(B)) is given by 



(N(B)) , (2.15) 



(N(B,s)) = -L £ NP(N,N b )B N s N » , (2.16) 

N,N b 

which can be considered as the average number of steps, made with fugacity B and stiffness 
s, by all possible SAWs that cross the fractal generator. Then, from ( 12.81) and ( 12.101) . it 
follows 

K = (N(B*, s)) . (2.17) 

The last formula enables us to calculate v via the MCRG method, that is, without 
calculating explicitly the coefficients P(N,Nf,). For a given fractal (with the scaling factor 
b) and the SAW stiffness s, we begin with determining the critical fugacity B*. To this 
end, we start the Monte Carlo (MC) simulation with an initial guess for the fugacity B in 
the region < Bq < 1. Here Bq can be interpreted as the probability of making the next 
step along the same direction from the vertex that the walker has reached, while sBq is 
the probability to make the next step by changing the step direction. We assume that the 
walker starts his path at one terminus (vertex) and tries to reach the other terminus of the 



TABLE I: The exact (3 < b < 9) and MCRG (11 < b < 201) results for the critical fugacities 
B* of the PF family of fractals. Each MCRG entry of the Table has been obtained by performing 
at least 10 5 MC simulations. The numbers in the brackets represent the MCRG errors concerning 
the last two digits, for instance, for b = 101 fractal the reading should be the following: B*(s = 
1) = 0.38815(05) = 0.38815 ± 0.00005. The values for b = oo are obtained by linear extrapolation 
of MCRG values. 



b 


B*(s = 1) 


B*(s = 0.9) 


B*(s = 0.7) 


B*(s = 0.5) 


B*(s = 0.3) 


B*(s = 0.1) 


3 


0.70711 


0.75595 


0.85923 


0.94815 


0.99212 


0.99990 


5 


0.59051 


0.63304 


0.73677 


0.86443 


0.97258 


0.99965 


7 


0.53352 


0.57132 


0.66544 


0.79312 


0.94262 


0.99924 


9 


0.50029 


0.53516 


0.62208 


0.74257 


0.90676 


0.99866 


11 


0.47863(23) 


0.51141(24) 


0.59352(26) 


0.70754(27) 


0.87251(13) 


0.99785(05) 


13 


0.46319(14) 


0.49449(20) 


0.57335(22) 


0.68289(24) 


0.84351(25) 


0.99679(07) 


15 


0.45191(13) 


0.48212(18) 


0.55857(20) 


0.66399(21) 


0.82073(22) 


0.99525(08) 


17 


0.44321(11) 


0.47307(17) 


0.54733(18) 


0.64984(19) 


0.80218(20) 


0.99315(09) 


21 


0.43065(10) 


0.45927(06) 


0.53091(07) 


0.62963(16) 


0.77512(17) 


0.98700(10) 


25 


0.42207(08) 


0.45005(12) 


0.51956(06) 


0.61549(14) 


0.75663(14) 


0.97767(11) 


31 


0.41323(06) 


0.44057(10) 


0.50811(11) 


0.60115(15) 


0.73802(12) 


0.96076(10) 


35 


0.40913(06) 


0.43611(09) 


0.50276(10) 


0.59441(11) 


0.72923(11) 


0.95010(10) 


41 


0.40420(04) 


0.43086(06) 


0.49664(04) 


0.58698(10) 


0.71913(10) 


0.93637(08) 


51 


0.39893(07) 


0.42498(07) 


0.48969(03) 


0.57825(15) 


0.70797(09) 


0.91977(07) 


61 


0.39524(06) 


0.42112(06) 


0.48493(07) 


0.57259(07) 


0.70032(08) 


0.90832(08) 


81 


0.39081(05) 


0.41638(05) 


0.47933(03) 


0.56545(07) 


0.69079(11) 


0.89413(05) 


101 


0.38815(05) 


0.41344(06) 


0.47601(11) 


0.56137(06) 


0.68541(06) 


0.88578(08) 


121 


0.38649(04) 


0.41167(07) 


0.47373(09) 


0.55862(05) 


0.68184(05) 


0.88007(04) 


151 


0.38479(02) 


0.40989(09) 


0.47163(04) 


0.55592(08) 


0.67825(05) 


0.87459(11) 


171 


0.38399(03) 


0.40894(06) 


0.47054(06) 


0.55464(04) 


0.67668(04) 


0.87217(07) 


201 


0.38316(06) 


0.40803(02) 


0.46951(11) 


0.55321(09) 


0.67480(04) 


0.86939(09) 


oo 


0.37915(40) 


0.40189(12) 


0.46217(15) 


0.54424(15) 


0.66287(22) 


0.85186(20) 
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TABLE II: The exact (3 < b < 9) and MCRG (11 < b < 201) results for the critical exponents 
v obtained in this work, for the PF family of fractals. Each MCRG entry of the Table has been 
obtained by performing at least 10 5 MC simulations. Numbers in the brackets correspond to the 
errors of the last two digits, determined by the MC simulation statistics (for example, for b = 11 



and s 


11 O TT7£1 |l f\TTn 

— U.y, WU IlcLVt! 




= 0.77239 ±0.00027). 






b 


u(s = l) 


u(s = 0.9) 


z/( s = 0.7) 


v(s — 5 s ) 


u(s = 0.3) 


v(s = 0.1) 


3 


0.79248 


0.81384 


fl 8793D 


n 944nn 


0.99061 


0.99988 


5 


0.78996 


0.79864 


82629 


88266 


0.96922 


0.99959 


7 


0.78111 


0.78666 


80342 


83910 


0.93233 


0.99906 


9 


0.77464 


0.77886 


791 08 


81 500 


0.88946 


0.99813 


11 


0.76959(27) 


0.77239(27) 


781 94f96") 


801 27(27") 


0.85478(15) 


0.99659(03) 


13 


0.76494(18) 


0.76923(25) 


D 77739f'95") 


fl 79950(95") 


0.83325(29) 


0.99415(11) 


15 


0.76232(17) 


0.76571(24) 


77307f24") 


78678(24") 


0.81825(26) 


0.99022(14) 


17 


0.75976(17) 


0.76226(23) 


76962f23") 


78205 ( 22") 


0.80885(24) 


0.98385(18) 


21 


0.75522(16) 


0.75754(10) 


7641 1 (09") 


77496(21 ") 


0.79717(21) 


0.96152(27) 


25 


0.75199(15) 


0.75406(07) 


D 76DD1 f09") 

U . 1 UUU1 1 UJ 1 


n 77nnn(i 9") 

U.I I UuU \JLcrj 


0.78993(20) 


0.92615(31) 


31 


0.74822(12) 


0.74954(20) 


7555914 9") 


76515(18") 


0.78315(18) 


0.87448(29) 


35 


0.74530(14) 


0.74773(19) 


75278(18") 


76276(47") 


0.77889(17) 


0.85207(25) 


41 


0.74332(08) 


0.74487(06) 


0.74966(08) 


0.75837(16) 


0.77439(16) 


0.83262(21) 


ol 


U. / oyoy(l i ) 


U. /414o(l 1 ) 


0.74614(07) 


0.75381(15) 


U. lb\)lU{Lo ) 


n oi inc. 1 1 1\ 
U.Ol 1 ZD(1 1 ) 


61 


0.73643(17) 


0.73930(16) 


0.74352(05) 


0.75051(14) 


0.76579(14) 


0.80992(15) 


81 


0.73314(15) 


0.73527(14) 


0.73935(06) 


0.74640(19) 


0.76050(12) 


0.80041(13) 


101 


0.73103(11) 


0.73280(04) 


0.73610(14) 


0.74313(12) 


0.75689(12) 


0.79482(12) 


121 


0.72865(13) 


0.73030(22) 


0.73459(28) 


0.74071(11) 


0.75405(11) 


0.79097(11) 


151 


0.72771(07) 


0.72840(06) 


0.73219(11) 


0.73848(10) 


0.75102(10) 


0.78618(08) 


171 


0.72631(12) 


0.72752(05) 


0.73094(10) 


0.73830(10) 


0.74850(09) 


0.78338(10) 


201 


0.72486(12) 


0.72633(05) 


0.72958(12) 


0.73583(04) 


0.74755(09) 


0.78077(10) 
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generator. In a case that the walker does not succeed to pass through the generator, the 
corresponding path is not taken into account. We repeat this MC simulation L times, for 
the same set B and s. Thus, we find how many times the walker has passed through the 
generator, and by dividing the corresponding number by L we get the value of the function 
( 12. 8p . denoted here by B'(B ,s). In this way we get the value of the sum 02. 8p without 
specifying the set P(N,Nb). Then, for a fixed s, the next values B n (n > 1), at which 
;he MC simulation should be performed, can be found by using the "homing" procedure 



251 ]. which can be closed at the stage when the difference B n — B n _i becomes less than the 
statistical uncertainty associated with B n ^\. Consequently, B* can be identified with the 
last value B n found in this way. Performing the MC simulation at the values B* and s, 
we can record all possible SAWs that traverse the fractal generator. Then, knowing such a 
set of walks, we can represent the average value of the length of a walk (that traverse the 
generator) via the corresponding average number of steps (N(B*,s)), and, accordingly, we 
can learn the value of the v through the formulas ( 12.171) and ( 12.9f) . In Tables [I] and HIl we 
present our MCRG results for B* and v respectively, for the chosen set of s values, for the 
PF fractal lattices with 11 < b < 201. 

To calculate the critical exponent 7 we need to find the singular behavior of the generating 
function G(x, s). The structure of the expression (12 ,4p shows that the asymptotic behavior of 
G(x, s), in the vicinity of the critical fugacity x c (s), depends on the corresponding behavior 
of the restricted partition functions (12.5l) -( !2T7j) . Assuming that the singular behavior of f]2 .41) 
is of the form ( 12.31) . it can be shown [18] that the critical exponent 7 should be given by 

7 = 2i«, (2.18) 
lnA^ 

where A 7 is the RG eigenvalue 

A 7 = a(5*,s), (2.19) 

of the polynomial a(B^ n ~ 1 \ s) defined by (12. 5p . with B* being the fixed point value of (12.81) . 
Therefore, it remains either to determine exactly an explicit expression for the polynomial 
a (f?*™ -1 ), s), or somehow to surpass this step and to evaluate only the single needed value 
a(B*,s). We have been able to determine the exact form of the requisite polynomial for 
the PF fractals with 3 < b < 9, while for b > 11 we have applied the MCRG to evaluate 
a{B*,s). 

In order to learn an explicit expression of the polynomial a(5^ r_1 - ) , s), we note that its 

10 



TABLE III: The exact (3 < b < 9) and MCRG (11 < b < 201) results for the critical exponents 
7 obtained in this work, for the PF family of fractals. Each MCRG entry of the Table has been 
obtained by performing at least 5 • 10 5 MC simulations. The numbers in the brackets represent 



the 


error bars related to the last two digits, for 


example for b = 


201 and s = 


0.9, we nave 


7 = 


2.216(13) = 2.216 ±0.013. 










b 


t(s — ll 


^(s — 91 


<-v(s — 71 


7(s = 0.5) 


7(s = 0.3) 


7(s = 0.1) 


3 


1 684n 

X . UOIU 


1 6796 


l.UUdu 


1.3368 


0.8189 


0.2524 


5 


1 7423 


1 7406 


1 79^fi 


1.6250 


1.1654 


0.3363 


7 


1 7614 

X . 1 UIt: 


1 7605 


1 7559 


1.7247 


1.4586 


0.4302 


9 


1 7807 


1 7795 


1 7748 

1. 1 1 t:0 


1.7596 


1.6335 


0.5340 


11 


1 8048 (321 


1 7987(311 

X. I UvJ l lOX I 


1 79D8(98l 

X . 1 JUG I iiO 1 


1.7753(25) 


1.7194(17) 


0.6498(07) 


13 
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form, due to the underlying self-similarity, should not depend on r, and, for this reason, in 
what follows we assume r — 1. Then, one can verify the following expression 

a(B, s)=J2 Q(N, N b )B N s N " , (2.20) 

N,N b 

where Q(N, N b ) is the number of all SAWs of iV steps, with N b bends, that start at any 
bond within the generator (r = 1) and leave it at a fixed exit. By enumeration of all relevant 
walks, the coefficients Q(N, N b ) can be evaluated exactly up to b = 9. For b = 3 fractal, the 
polynomial (I2.20p is of the form 

a(B,s) = 1 + B + IBs + B 2 + 2B 2 s + 2B 2 s 2 + 2B 3 s 2 
+ 2B 3 s 3 + AB A s 3 + 4BV + 25V + 2B 6 s 5 , 

(2.21) 

for b = 5 it is given in the Appendix, while for 6 = 7 and 9 they are given in the supple- 



mentary EPAPS Document [27|]. Using this information, together with (12.191) . (12.181) . and 
previously found B* and A„, we have obtained the desired exact values of 7 (see Table IITTI) . 

For a sequence of b > 9, the exact determination of the polynomial (12.201) . that is, knowl- 
edge of the coefficients Q(N,N b ), can be hardly reached using the present-day computers. 
However, to calculate A 7 one does not need a complete knowledge of polynomial a(B, s), but 
only its values at the fixed point (see Eq. (I2.19P ). However, the polynomial that appears in 
(12. 5p can be conceived as grand partition function of an appropriate ensemble, and conse- 
quently, within the MCRG method, the requisite value of the polynomial can be determined 
directly jis| . Owing to the fact that we can obtain A 7 = a(B* } s) through the MC simula- 
tions, and, knowing X u from the preceding calculation of u, we can apply (I2.18P to calculate 
7. In Table ITO1 we present our MCRG results of 7 for 11 < b < 201, for the chosen set of 
stiffness parameter values (s = 1, 0.9, 0.7, 0.5, 0.3 and 0.1). 



III. DISCUSSION AND SUMMARY 



We have studied critical properties of semiflexible polymer chains on the infinite family 
of the PF fractals whose each member has the fractal dimension df equal to the Euclidean 
value 2. In particular, we have calculated the critical exponents v and 7 via an exact RG 
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FIG. 3: (Color online) Plot of the end-to-end distance critical exponent v of semiflexible polymer 
chain on PF fractals: (a) as a function of stiffness parameter s, for various fractal scaling parameter 
b, and (b) as a function of 1/6, for various values of stiffness parameter s (the horizontal dashed 
line represents the Euclidean value v = 3/4 for flexible polymers, s = 1, whereas thin solid lines 
serve only as guides to the eye). 

(for 3 < b < 9) and via the MCRG approach (up to b = 201). Specific results for the critical 
exponents have been presented in Tables [III and IHII 

In order to analyze the obtained results, in Fig. [3](a) we have plotted v as a function of 
stiffness parameter s, for several values of fractal scaling parameter b. One can see that for 
each b, exponent v monotonicalfy decreases from the value v = 1, for s = 0, corresponding 
;o rigid rod, to the value vsAwip), for s = 1, corresponding to the flexible polymer chain 
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281 ] . This, indeed, implies that for finite b, the mean end-to-end distance for semiflexible 
polymers increases with its rigidity, and is always between its values for the flexible chain 
and the rigid rod. In the same figure one can also observe that when b increases the curves 
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v(s) become increasingly sharper so that their limit looks to be v = 1, at s = 0, whereas 
v ~ const., for < s < 1. This observation may imply that for very large b (beyond b = 201) 
the critical exponent v becomes independent of s. Here, one should note that it is believed 
that critical exponent v is universal for semiflexible SAWs on Euclidean lattices, that is, it 
does not depend on s This expectation is based on universality arguments, and it was 
exactly demonstrated for directed semiflexible SAWs [jjj]. The same conclusion was also 
exactly derived for semiflexible SAWs on the Havlin-Ben-Avraham and 3-simplex fractal 
lattices [ijj] . However, as it was pointed out in [ijj] , in contrast to the case of homogeneous 
lattices, where rigidity only increases the persistence length, but does not affect neither 
the scaling law governing the critical behavior of the the mean end-to-end distance, nor 
the value of the critical exponent v of SAWs, one might expect that presence of disorder 
in nonhomogeneous lattices, combined with the stiffness, in some cases can constrain the 
persistence length, and consequently induces dependence of ^ on s. This was explicitly 
confirmed in the same paper, by exact calculation of the critical exponent v for branching 
Koch curve, which turned out to be continuously decreasing function of s, similar to functions 
depicted in Fig. [3]^a). Apparently, the established dependence of v on s for PF fractals with 
smaller values of b shows that considerable lattice disorder affects significantly the values of 
v, while the dependence of v on s gradually disappears for PF fractals with smaller disorder 
(appearing for larger b). These facts confirm the assumption [15( that lattice disorder, 
combined with the polymer stiffness, has a predominant impact on the critical behavior of 
semiflexible polymers. 

In Fig.[3](b) data for v as a function of 1/6 are depicted, for various values of s. It appears 
that for each considered value of stiffness s in the range < s < 1, the critical exponent v 
is monotonic function of the scaling parameter b, in the region of b studied. It can be also 
seen that for large fixed b, the differences between the values of v(s) (for various s) decrease 
when b increases, which brings us to the question of the behavior of v in the fractal-to- 
Euclidean crossover, when b — > oo. Concept of the fractal-to-Euclidean crossover is often 
used in order to study if and how various physical properties change when inhomogeneous 
lattices approach homogeneous (translationally invariant Euclidean) lattice. By tuning some 
conveniently chosen parameter of the fractal lattice (such as scaling parameter b in the case 
of PF fractals), properties of the corresponding Euclidean lattice (square lattice in this case) 
can be gradually approached. Studies of the flexible SAW models on Sierpinski gasket family 
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of fractals 
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26 



29l . |30j, |3l|, l32j, |33j, as well as on PF fractals 18), revealed that crossover 



behavior of critical exponents can be rather subtle in the sense that not all critical exponents 
tend to their Euclidean values, and even when they do so it can be accomplished in quite 
unexpected manner. For instance, according to the finite-size scaling arguments, when 
b — > oo exponent v of flexible polymers (s = 1) on PF fractals approaches the Euclidean 
value 3/4 from below 18j|, which together with the fact that v is monotonically decreasing 
function for b up to 201, means that for some value of b larger than 201 there should exist a 
minimum. On the other hand, for s = exponent v is equal to 1 for each b. For < s < 1 
apparent trend of the curves presented in Fig. [3]^b) suggests that limiting value of z/, when 
b — > oo, does not depend on particular value of s, and following the behavior of v for flexible 
polymers, it should be equal to the Euclidean value 3/4. However, we would not like to 
draw here such a definite conclusion without additional investigations. 

Continuing the comparison of the criticality of flexible and semiflexible SAWs on the PF 
fractals, in Fig. [4] we have depicted the data from Table [I] for the critical fixed points B* . On 
the left-hand side of this figure one can notice that B*, which is equal to the reciprocal of the 
connectivity constant /i, is monotonically decreasing function of s, for each b considered. This 
has been expected, since /i has the physical meaning of the average number of steps available 
to the walker having already completed a large number of steps, so that larger flexibility of 
the polymer chain implies larger /i, and consequently fi(s = 0) < /i(0 < s < 1) < fi(s = 1). 
In Fig. H](b) one can also observe that for fixed s, the fixed point B* decreases with b. 
Furthermore, B* becomes almost linear function of 1/6 for large b, which allows us to 
estimate the limiting values of B* for b — > oo. The obtained asymptotic values are given at 
the end of Table H The value B*(s = 1) = 0.37915 ± 0.00040, should be compared with the 
Euclidean value 0.3790523(3) for the square lattice, obtained in 34|. As one can see, the 
agreement is very good, and we may say that for flexible polymers the values B*(b), in the 
limit b — > oo, converge to the Euclidean value. Similarly, we expect that, for given s < 1, the 
values B*(b) of semiflexible polymers also converge to the corresponding d = 2 Euclidean 
values (which are functions of s). In Fig. [5](a) we have plotted the lines obtained by linear 
fitting of the large b data for B*(s) for s = 0.1, 0.3, 0.5, 0.7, 0.9 and s — 1. In the present 
situation, estimated limiting values of the connectivity constants fi(s) = 1/B*[s, b — > oo) are 
depicted in Fig.^b), as function of s. It seems that n(s) is linear function of s, implying that 
connectivity constant for semiflexible SAWs on square lattice could be a linear function of 



15 




FIG. 4: (Color online) Data for the fixed point values B* (the reciprocal connectivity) of semiflex- 
ible polymer chain plotted as: (a) function of stiffness s, for various values of the fractal parameter 
6, and (b) function of 1/6, for various values of s. 

the stiffness parameter s. Such expectation is also in accord with the exact results obtained 
for directed semiflexible SAWs [14|. 

To make our analysis of semiflexible SAWs on PF fractals complete, in Fig. [6] we present 
the data found for the critical exponent 7. As it was explained in Sec. [TTJ exponent 7 is given 
by (12.181) : 7 = 2hi(A 7 /6)/lnA l/ , where eigenvalues X u and A 7 (given by (12.101) and (12.191) . 
respectively) are evaluated at the fixed point B = B* of the RG equation (I2.6p . For 6 = 3 
the fixed point B* is given by (12.131) . and since dependence of X u and A 7 on B is known 
by exact means, the curve 7(6 = 3, s) in Fig. [6](a) was plotted according to the closed-form 
exact formula. For 6 = 5, 7 and 9, RG equation ( 12. 6ft was also found explicitly, but its 
fixed point can be calculated only numerically in these cases. Nevertheless, this task can 
be done for large number of s values, and putting fixed points B*(s) calculated accordingly, 
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FIG. 5: (Color online) (a) Linear fitting of the data (full squares) for the fixed point B* (see Table|I]) 
as function of 1/6 for large values of b, and various s. Circles correspond to the extrapolated limiting 
values of B*(s, b — ► oo), and their precise values are given in the last row of Table [IJ (b) Plot of the 
connectivity constant fi(s) = 1/B*(s,b — ► oo). Dashed line is the linear fit of the data presented 
by circles, which correspond to /x(s) for s = 0, 0.1, 0.3, 0.5, 0.7, 0.9 and s = 1. 

into the exact expressions found for X u and A 7 , one obtains the corresponding values for 7, 
and, consequently, curves for b = 5, 7 and 9 in Fig. E(a). For larger values of b, depicted 7 
curves were obtained by interpolating the data found by MCRG approach for s = 0.1, 0.3, 
0.5, 0.7, 0.9 and s = 1 (Table HTT1) . and generalizing the fact 7(6, s = 0) = 0, exactly found 
for smaller values of b, to all b values. One can see that for each b, the critical exponent 7 is 
monotonically increasing function of the stiffness parameter s. The dependence of 7 on s is 
in accord with the discussed non-universality of u, and also with the results obtained for 7 
of SAWs on the branching Koch curve and Havlin-Ben-Avraham fractal [15J. However, one 
should note here that while for the Koch curve both v and 7 depend on s, in the case of 
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FIG. 6: (Color online) Data for the SAW critical exponent 7 presented as: (a) function of the 
stiffness parameter s, for various values of the fractal enumerator b, and (b) function of 1/6, 
for s = 1, 0.9, 0.7, 0.5, 0.3, and 0.1 (the horizontal dashed line represents the two-dimensional 
Euclidean value 7 = 43/32 for flexible polymers s = 1). 

Tavlin-Ben-Avraham fractal only 7 is non-universal [y = 1 for all values of s). Besides, in 
l^ l it was shown that neither v nor 7 depend on s for SAWs on the 3-simpiex fractal. The 
observed different behavior of exponents v and 7 on various fractals is an intriguing fact 
and imposes the question of the universality of 7 for semiflexible SAWs on homogeneous 
lattices. One might try to draw a helpful conclusion by looking at the large b behavior of 
the functions 7(5), plotted in Fig. [6](a). One can see that as b grows the curve j(s) becomes 
sharper, so that for b = 201 it is almost constant in the large part of the region < s < 1, 
whereas in the vicinity of s — it rapidly drops to the value 7 = 0, at s = 0. Therefore, 
it may be concluded that for b ^> 1 and < s < 1 exponent 7 becomes independent of 
s. Furthermore, in Fig. E^b), we perceive that for each studied s, the critical exponents 7 
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monotonically increase with b, and for b = 201 acquire almost the same value 7 ~ 2.2. These 
observations may imply that 7 for semiflexible SAWs on homogeneous lattices is universal. 
However, it is known that critical exponent 7 for flexible polymers (s = 1) on the two- 
dimensional Euclidean lattices is equal to 7 = 43/32, which is far from the apparent limiting 
value 2.2 (suggested by the plots in Fig.H^b), when 1/6 —> 0), implying that 7 for PF fractals 
does not tend to its Euclidean value for large b. This may seem odd, but it fits quite well 
into the peculiar picture which have emerged during the last two decades for the fractal-to- 
Euclidean crossover behavior of the exponent 7 of flexible SAWs on PF 18| and on Sierpinski 
gasket (SG) family of fractals [20j, |32|, |33j , as well as for some models of directed SAWs on SG 
fractals [31I - For instance, using finite-size scaling arguments, Dhar [2^] concluded that 7 
for SAWs on SG fractals at the fractal-to-Euclidean crossover approaches the non-Euclidean 
value 133/32. In a similar manner, for PF fractals it was also demonstrated [18( that in 
the limit b — > 00 exponent 7 tends to 103/32, which is again the non-Euclidean value. In 
addition, numerical analysis of the large set of exact values of 7 obtained for the piece-wise 
directed SAWs on SG fractals, as well as an exact asymptotic analysis [3JJ, also showed that 
the limiting value of 7 differs from the corresponding Euclidean value. In all these cases the 
established crossover behavior could not have been predicted only on the basis of 7 values 
obtained for relatively small b (up to b = 201, for instance). On these grounds we may 
infer that in the crossover region, when b — > 00, critical exponent 7 does not depend on the 
stiffness s, and approaches the non- Euclidean value. 

In conclusion, we may say that family of plane-filling fractals proved to be useful for 
investigation of the effects of the rigidity on the criticality of SAWs on nonhomogeneous 
lattices. It is amenable to applying exact and Monte Carlo renormalization group study, 
which we performed on large number of its members. The obtained results show that 
the critical behavior of semiflexible SAWs is not universal, in a sense that critical end- 
to-end exponents is, as well as entropic exponents 7, continuously vary with the stiffness 
parameter s. Such non-universality does not occur on regular lattices, but it was found for 
SAWs on the branching Koch curve, suggesting that polymer behavior in realistic disordered 
environment might be more affected by its stiffness than it was expected. Apart from the 
stiffness parameter s, critical exponents also depend on the fractal parameter b, but the 
trend of functions is(s) and 7(5) is similar for different b values. This similarity becomes 
more pronounced as b grows and approaches the fractal-to-Euclidean crossover region (b — ► 



19 



oo). Assuming that critical exponents on regular lattices do not depend on s, it would be 
challenging to reveal what exactly happens with the exponents in the limit b — > oo, which 
we would like to investigate in the future. 
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APPENDIX 

Here we give the coefficients Q(N, Nj,) of the RG equation (I2.20p for b = 5 PF fractal: 
Q(0, 0) = 1, Q(l, 0) = 1, Q(l, 1) = 2, Q(2, 0) = 1, Q(2, 1) = 2, Q(2, 2) = 2, Q(3, 0) = 1, 
<5(3, 1) = 4, Q(3, 2) = 6, Q(3, 3) = 4, Q(4, 0) = 1, Q(4, 1) = 4, Q(4, 2) = 10, Q(4, 3) = 12, 
Q(4,4) = 6, Q(5,2) = 6, Q(5,3) = 18, Q(5,4) = 18, Q(5,5) = 6, Q(6,2) = 2, Q(6,3) = 20, 
Q(6,4) = 38, Q(6,5) = 32, Q(6, 6) = 12, Q(7, 2) = 2, Q(7, 3) = 8, Q(7,4) = 28, Q(7, 5) = 
34, Q(7,6) = 32, Q(7,7) = 10, Q(8,3) = 4, Q(8,4) = 16, Q(8,5) = 52, Q(8,6) = 62, 
Q(8,7) = 48, Q(8,8) = 14, Q(9,4) = 10, Q(9, 5) = 34, Q(9, 6) = 70, Q(9, 7) = 54, 
Q(9,8) = 34, Q(9,9) = 6, Q(10,5) = 18, Q(10,6) = 38, Q(10,7) = 78, Q(10,8) = 68, 
Q(10,9) = 44, Q(10,10) = 4, Q(ll,5) = 2, Q(ll,6) = 14, Q(U,7) = 44, Q(ll,8) = 70, 
Q(ll,9) = 52, Q(ll,10) = 36, Q(ll,ll) = 6, Q(12, 7) = 18, (5(12,8) = 30, (5(12,9) = 56, 
<5(12, 10) = 42, Q(12, 11) = 34, Q(12, 12) = 6, Q(13, 8) = 4, Q(13, 9) = 20, Q(13, 10) = 28, 
Q(13, 11) = 30, g(13, 12) = 22, Q(13, 13) = 2, Q(U, 11) = 22, Q(U, 12) = 22, Q(U, 13) = 
12, Q(15,12) = 4, Q(15,13) = 8. 
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